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ABSTRACT 

A  method  for  the  calculation  of  the  reaction 
matrix  of  the  Brueckner  theory  is  described.  The  two- 
nucleon  interaction  is  separated  into  two  parts.  The 
separation  distance  is  so  defined  that  the  energy  is 
the  eigenvalue  of  the  effective  Hamiltonian  which 
consists  of  the  long  range  part  of  the  interaction  and 
the  kinetic  energy.  The  short  range  interaction  is 
treated  as  a  perturbation  which  modifies  the  wave  function. 
If  the  outer  potential  is  weak,  a  plane  wave  approximation 
is  used  for  the  unperturbed  wave  function  and  the 
diagonal  elements  of  the  reaction  matrix  are  given 
to  first  order  by  the  first  Born  approximation  of  the 
long  range  potential.  The  separation  distance  is 
calculated  self-consistently,,  the  iteration  being 
simplified  by  the  used  of  the  "reference  spectrum" 
approximation.  Finally  there  is  a  discussion  of  the 
singularities  in  the  integral  equations  of  the 
Brueckner  theory  and  the  divergence  of  the  expansion 
of  the  reaction  matrix.  Attempts  to  modify  the  theory 
so  as  to  eliminate  these  difficulties  are  described. 
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Chapter  I.  Introduction 


The  first  difficulty  encountered  in  attempting 
a  theoretical  study  of  the  properties  of  nuclear  matter 
is  the  representation  of  the  various  interactions 
present  in  the  system.  A  complete  description  would  require 
a  formulation  in  terms  of  second  quantized  fields  containing 
terms  involving  the  co-ordinates  of  two,  three,  and 
arbitrarily  many  nucleons  and  terms  describing  the  interaction 
of  nucleons  with  a  meson  field.  Such  a  system  is  very 
complicated  and  it  is  obviously  necessary  to  construct 
a  model  which  is  mathematically  tractable.  The 
traditional  (and  most  successful)  approach  has  been  to 
eliminate  the  meson  field  and  the  many-body  interactions, 
leaving  an  effective  Hamiltonian  consisting  of  just  two- 
particle  interactions  plus  the  kinetic  energy.  Another 
problem  then  arises,  viz.,  expressing  the  two-particle 
interaction  as  a  potential  model.  It  is  well  known  that 
a  potential  model  can  be  justified  only  for  the  description 
of  low-energy  phenomena  in  which  the  separation  of  the 
two  nucleons  is  at  least  one  half  the  pion  Compton  wave¬ 
length,  i.e.,  for  the  long  and  medium  ranges  of  the 
interaction  as  described  by  Taketani  et  al ^ ^ .  The  long 
range  potential  is  well  known.  It  is  the  asymptotic  part 
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’ 
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of  the  one-pion  exchange  potential  of  meson  theory.  In 

the  medium  range  one  and  two-meson  exchanges  contribute 

and  the  potential  is  not  known  accurately  but  it  is  possible 

to  represent  the  interaction  by  a  phenomenological  static 

( 2 ) 

potential  .  At  short  range ,  exchanges  of  many  resonances 
and  pion  resonances  contribute  and  relativistic  and 
recoil  effects  are  substantial.  In  this  region  the  use 
of  a  potential  model  may  not  be  valid  even  as  an  approxi¬ 
mation.  The  complete  potential  must  produce  the  correct 
two  nucleon  scattering  phase  shifts  but  this  does  not 

determine  the  potential  unless  the  potential  is  assumed 

(q  4) 

to  be  static v  .  Even  then  the  model  might  not  be 

adequate  for  calculating  the  properties  of  nuclear 

( 5) 

matter v ^ ' . 


The  two-nucleon  interaction  thus  consists  of  a 

long  range  part  which  is  fairly  well  known  and  can  be 

expressed  by  a  static  potential  and  a  short  range  part 

which  in  practice  can  be  represented  by  a  potential 

model  or  by  a  boundary  condition  on  the  two-nucleon  wave 

function.  In  the  innermost  region  the  interaction 

appears  to  be  very  strongly  repulsive.  This  is 

usually  represented  by  introducing  the  so-called  "hard 

core"  and  making  the  potential  infinity  positive  at 

-13 

distances  less  than  about  0.4  x  10 


cm.  Outside  the 
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hard  core  radius  the  potential  is  attractive.  In 
addition  a  tensor  force  is  required  to  explain  the 
quadrupole  moment  of  the  deuteron  and  spin-orbit  and 
quadratic  spin-orbit  interactions  are  probably  necessary 
to  fit  high  and  low  energy  scattering  data. 

Electron  scattering  experiments  have 

shown  that  the  mean  spacing  of  nucleons  in  nuclear 

-13 

matter  is  about  1.8  x  10  J  cm  which  is  large  compared 

with  the  radius  of  the  hard  core.  Thus  the  effect  of  the 

hard  core  may  be  regarded  as  weak  and  the  ideal  fermi 
gas  state*  which  is  the  starting  point  for  nearly  all 
current  theories*  may  be  a  good  approximation  to 
the  actual  state  of  the  system*  depending  on  the 
strength  of  the  rest  of  the  potential.  The  proposition 

that  the  system  resembles  a  gas  more  than  a  liquid  is 

supported  experimentally  by  the  fact  that  low  energy 
nucleons  incident  on  nuclei  have  mean  free  paths  which 
are  large  compared  with  the  mean  separation  of  nucleons. 

Two  important  methods  have  been  used  to 
calculate  the  properties  of  nuclear  matter.  One  method 
is  based  on  the  variational  principle  and  the  use  of 
trial  wave  functions.  Use  of  the  Slater  determinental 
wave  function  in  which  the  single  particle  states  are 
plane  waves  with  momenta  less  than  or  equal  to  the 
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fermi  momentum  (i.e.  the  ideal  fermi  gas  state)  as  a 
trial  wave  function  gives  an  upper  limit  of  (+  °°)  for 
the  energy  per  particle.  One  way  to  avoid  this  spurious 
result  caused  by  the  hard  core  is  to  impose  a  boundary 
condition  on  the  trial  wave  function,  in  particular  to 
multiply  it  by  a  correlation  function  which  vanishes 
whenever  two  particles  approach  within  the  hard  core 
radius.  This  method  was  proposed  by  Jastrow  (6).  The 
procedure  is  quite  laborious  and  no  satisfactory  results 
have  been  obtained  from  first  order  calculations  (7). 

The  most  successful  approach  to  understanding 
the  theory  of  nuclear  structure  is  the  Brueckner  theory 
which  is  a  development  of  the  perturbation  method. 

The  basic  idea  is  that  the  energy  may  be  obtained  by 
replacing  the  potential  by  an  operator  G  which  is 
generated  by  an  integral  equation  (equation  (3*29)  )• 

This  operator  contains  the  effects  of  two-body  collisions 
in  all  orders  of  perturbation  theory  while  ignoring 
many-body  effects.  The  operator  G  differs  from  the 
reaction  matrix  of  scattering  theory  only  in  the  presence 
of  the  exclusion  principle  and  a  different  energy 
denominator  and  because  of  this,  it  is  usually  referred  to 
as  the  Brueckner  reaction  matrix  or  just  the  reaction 
matrix.  Now  the  reaction  matrix  of  scattering  theory 
is  finite  even  when  the  potential  has  a  hard  core  and 
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each  term  in  its  perturbation  series  is  meaningless.  Hence 
it  appears  reasonable  to  suppose  that  the  Brueckner  reaction 
matrix  is  also  well  defined  and  finite  for  singular 
potentials.  The  convergence  properties  of  the  perturbation 
series  are  of  course  unknown  and  the  above  assumption  is 
not  justified  mathematically.  It  is  customary  to  ignore 
this  difficulty  and  proceed  blithely  with  the  calculations. 
It  is  now  known,,  however,  that  the  Brueckner  matrix  is 
not  always  convergent.  This  difficulty  is  the  subject 
of  chapter  5* 

The  actual  calculation  of  the  Brueckner  matrix 
is  a  very  difficult  problem  because  of  the  self-consistency 
requirements  on  the  solution.  The  energy  denominator  in 
the  integral  equation  for  G  contains  a  potential 
energy  contribution  which  depends  on  G  and  it  is  the 
coupling  of  the  two  equations  which  complicates  the 
problem.  Numerical  solutions  can,  in  theory,  be  obtained 
by  iteration.  In  practice,  approximations  are  required. 

A  comprehensive  discussion  is  given  by  Brueckner  and 
Gamine  1  (8)  who  used  a  realistic  potential  and  obtained 
good  results  for  the  binding  energy  per  nucleon  and  the 
density  of  nuclear  matter. 


( 
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However  the  numerical  calculations  required  for 
an  iterative  solution  of  the  Brueckner  equations  are 
complicated  and  not  very  transparent.  It  is  thus 
desired  to  simplify  the  method  of  applying  the  Brueckner 
theory.  Important  progress  in  this  direction  was  made 
by  Moszkowski  and  Scott  (9).  They  separated  the 
potential  into  a  short  range  part  and  a  long  range  part, 
the  point  of  separation  being  such  that  the  short  range 
potential  gives  no  phase  shift  for  free  particle 
scattering.  The  short  range  potential  is  then  represented 
by  a  short  range  reaction  matrix  which  vanishes  for  free 
nucleons  (with  interactions  between  the  two  particles  but 
no  interaction  with  the  rest  of  the  system)  and  is 
small  for  nuclear  matter  since  the  two  particle  wave 
function  in  nuclear  matter  is  closely  approximated  by  the 
free  particle  wave  function  at  small  distances.  For 
larger  distances  the  wave  function  can  be  approximated 
by  the  unperturbed  wave  function  and  the  contribution  of 
the  long  range  potential  can  be  treated  by  Born 
approximation.  Thus  the  reaction  matrix  consists  of  the 
first  Born  approximation  for  the  long  range  potential  and 
the  difference  between  the  short  range  reaction  matrices 
for  nuclear  matter  and  for  free  nucleons,  plus  some  small 
corrections.  This  method  is  much  simpler  than  a  direct 
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solution  of  the  Brueckner  equations  but  Kohler  (10)  has 
shown  that  the  Moszkowski  -  Scott  calculation  of  the 
short  range  contribution  may  not  be  very  accurate. 

Another  method  of  simplifying  the  calculation 

of  the  reaction  matrix  of  the  Brueckner  theory  is  the 

reference  spectrum  method  introduced  by  Bethe  et  al  (ll). 

They  replace  the  single  particle  energy  spectrum  for  the 

intermediate  states  by  a  reference  spectrum  of  the 

2 

form  A  +  Bk  where  k  is  the  single  particle 
momentum.  Then  the  energy  denominator  will  be  a  simple 
quadratic  function  of  the  relative  momentum  of  the 
excited  state.  If  the  exclusion  operator  is  also  dropped 
the  integral  equation  can  be  reduced  to  a  second  order 
differential  equation  which  can  be  handled  quite  easily. 
The  reference  spectrum  parameters ,  A  and  B,  are  chosen 
to  fit  the  actual  particle  energies  with  the  reference 
spectrum  as  accurately  as  possible  for  k  >  kp.  The 

p 

reference  matrix  G  is  then  used  to  calculate  the 
Brueckner  matrix  G.  A  simple  quadrature  is  found  to  be 

■p 

sufficient.  The  difference  G  -  G  is  due  mainly  to  the 
effect  of  the  exclusion.  The  spectral  correction,  which 
arises  from  the  approximation  to  the  energy  denominator, 
is  quite  small  and  shows  that  a  detailed  knowledge  of 
the  form  of  the  single  particle  energy  spectrum  is  not 
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important  for  the  calculation  of  the  reaction  matrix. 

The  actual  energy  spectrum  must  of  course  be  calculated 
from  the  reaction  matrix. 

It  is  our  intention  here  to  develop  a  method 
of  calculating  the  Brueckner  reaction  matrix  which 
minimises  the  effect  of  the  purely  phenomenological 
part  of  the  potential,  i.e.,  the  innermost  region. 

The  starting  point  is  the  Hamiltonian  used  by  Brueckner, 
i.e.,  containing  two-body  interactions  only.  We 
separate  the  two-body  interaction  into  a  short  range 
part  and  a  long  range  part  and  treat  the  contributions 
to  the  Hamiltonian  of  the  short  range  interactions  as 
a  perturbation.  As  Moszkowski  and  Scott  (9)  have  shown 
the  Brueckner  matrix  is  given  to  first  order  by  the 
sum  of  two  terms  which  are  obtained  from  the  short  range 
potential  and  long  range  potential  respectively.  Cross 
terms  are  included  in  the  higher  orders  (c.f.  equation 
(3.34)  ).  We  choose  to  define  the  separation  point  by 
requiring  that  the  first  order  contribution  of  the 
short  range  potential  to  the  reaction  matrix  be  zero. 
This  differs  from  the  Moszkowski -Scott  definition  of 
the  separation  distance.  With  their  definition  the 
short  range  potential  contributes  to  the  first  order 
energy  shift,  the  contribution  being  small  if  the  free 
particle  wave  function  is  approximately  the  same  as  the 
nuclear  wave  function  at  small  distances.  On  the  other 
hand,  with  our  separation  the  first  order  expression 
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for  the  reaction  matrix  is  just  the  first  Born  approximation 
for  the  long  range  part  of  the  potential  and  the  energy 
of  the  system  can  be  obtained  from  the  effective 
Hamiltonian  which  contains  the  long  range  interaction 
and  the  kinetic  energy. 

If  only  two-body  correlations  are  considered,, 
the  problem  can  be  reduced  to  the  problem  of  the 
scattering  of  two  nucleons  in  the  field  of  the  average 
potential  produced  by  the  other  nucleons.  Since  the 
effective  Hamiltonian  contains  only  the  long  range 
potential,  the  equations  for  its  eigenfunctions  (the 
unperturbed  single  particle  wave  functions)  can  be 
solved  by  the  Hartree  -  Fock  method  and  the  solutions 
will  determine  the  effective  Hamiltonian  and  the 
exact  wave  function  of  the  system.  Such  a  solution  must 
be  self-consistent  of  course  and  the  procedure  would 
be  tedious.  The  problem  can  be  greatly  simplified  if 
we  assume  that  the  unperturbed  wave  functions  can  be 
approximated  by  plane  waves.  Then  the  first  order 
calculation  of  the  separation  distance  can  be  performed 
using  only  the  known  two-nucleon  scattering  phase  shifts 
and  the  form  of  the  long  range  part  of  the  potential, 
the  method  being  similar  to  that  used  by  Moszkowski  and 
Scott  (9).  From  this  result  we  can  calculate  an 
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approximate  value  for  the  single  particle  potential 
energy  and  this  in  turn  enables  us  to  obtain  a  more 
accurate  expression  for  the  separation  distance  if  the 
short  range  potential  is  known.  This  iteration  can  be 
continued  until  a  self-consistent  result  is  obtained, 
i.e.  until  the  difference  between  two  successive  results 
becomes  negligible. 

The  details  of  this  method  of  solving  the 
Brueckner  equation  are  contained  in  the  next  three 
chapters.  Chapter  2  contains  a  brief  outline  of  the 
two-body  scattering  problem  which  is  formulated  in  such  a 
way  that  the  equations  are  analogous  to  those  obtained 
in  an  approximate  solution  to  the  many-body  problem. 

The  equations  obtained  are  useful  in  suggesting  a  simple 
method  of  obtaining  the  separation  distance  for  the 
many-body  problem.  Chapter  3  contains  the  formal  develop¬ 
ment  of  our  method  of  attacking  the  many-body  problem 
and  shows  that  the  reaction  matrix  (and  therefore 
quantities  of  interest  such  as  binding  energy  and 
single  particle  energy)  can  to  a  very  good  approximation 
be  obtained  from  a  knowledge  of  the  long  range  potential 
and  the  separation  distance.  The  rearrangement  energy 
and  the  saturation  condition  are  discussed  in  terms  of 
the  derivatives  of  the  separation  distance.  Chapter  4  is 
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a  discussion  of  methods  of  calculating  the  separation 
distance  with  particular  emphasis  on  a  method  which 
utilizes  the  reference  spectrum  approximation  described 
above.  The  equations  then  reduce  to  those  of  chapter  2 
and  can  be  easily  solved.  First  order  corrections  to 
this  reference  spectrum  result  are  discussed.  The 
basic  ideas  contained  in  chapters  3  and  4  have 
previously  been  included  in  a  research  paper  (12). 

A  more  detailed  account  is  given  here. 

As  previously  mentioned  the  convergence  of 
the  series  expansion  of  the  Brueckner  reaction  matrix 
is  usually  assumed.  It  is  fairly  easy  to  show  however 
that  there  are  regions  in  which  the  expansion  is  divergent. 
Chapter  5  is  a  discussion  of  this  difficulty  and  of 
attempts  to  eliminate  the  divergence  by  applying  the 
Bardeen  -  Cooper-Schrief fer  theory  of  superconductivity. 
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Chapter  2.  The  Two  Body  Problem 

Before  any  equations  are  written  it  should  be 

noted  that  all  energy  quantities  are  defined  as  the 

_2 

actual  energy  multiplied  by  M  'fi  where  M  is  the 
mass  of  a  particle  of  the  system  under  consideration. 

Also  the  fermi  (l  P  =  10  3  cm)  is  used  as  a  unit  of 

-2 

length  so  that  energies  are  given  in  F  .  If  M  is 
the  nucleon  mass.,  one  may  translate  back  to  familiar 
units  by  using  the  relation 

1  F-2  =  41.467  Mev  . 

Consider  now  the  Hamiltonian  for  a  system  of 
two  particles  (of  equal  mass)  interacting  through 
the  central  potential  v(r), 

H  =  T  +  v (r )  ,  (2.1) 

where  T  is  the  kinetic  energy  operator  and  v(r) 
is  assumed  to  have  the  form  of  the  two-nucleon  potential, 
i.e.,  a  hard  core  of  radius  c  with  an  attractive 
interaction  for  r  >  c.  The  Hamiltonian  is  separated 
into  two  parts 

H  ~  H  +  H*  (2.2) 

o 

H  =  T  +  v(r)  0 (r  -  £ ) 
o 


where 


(2.3) 
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and  H’  =v(r)0(|-r)  (2.4) 

where  0(x)  is  the  step  function  equal  to  unity  for 
x  >  0  and  zero  for  x  <  o.  The  parameter  |  will 
be  specified  later. 

Let  |  if/  >  be  the  ground  state  eigenket  of 
the  total  Hamiltonian  corresponding  to  energy 
eigenvalue  e 


H  I  ^  >  =  e  |  ^  >  , 


(2.5) 


and  let  the  unperturbed  Hamiltonian  Hq  determine  the 
complete  orthonormal  set  of  eigenkets  |  b^>  correspond¬ 
ing  to  eigenvalues  e^ 


With  the  normalization 
|  tfj  >  in  terms  of  the  set 


=  e  |  b  >  .  (2.6) 

n  n 

<b  |  =  1,  the  expansion  of 

b  >  yields 

I  n  cX 


|  b  ><b  |  HT  |  ip  > 

ip  >  =  |  b  >  +  2  - 

n*0  e  -  en 


(2.7) 


and 


=  e  +  <b  |  H’  |  f  > 
o  o 


e 


(2.8) 
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The  separation  distance  £  is  now  chosen  to  be 
such  that 


<b  h»  \f>  =  0 
o  1  r 


(2.9) 


and 


e 


(2.10) 


It  is  assumed  that  the  potential  is  such  that  a  finite 
£  exists.  As  defined  by  equation  (2.9)*  £  will  be 
a  function  of  the  energy  eQ  or  of  the  magnitude  kQ 
of  the  relative  momentum  of  the  two  particles  if  the 
|  b^  >  are  momentum  states. 


If  the  long  range  part  of  the  potential  is 
weak,,  the  eigenvectors  of  Hq  may  be  approximated 
by  the  free  particle  wave  vectors  which  are  represented 
in  co-ordinate  space  by  plane  waves,  i.e. 


and 


<  r 


b  >  ai 
n 


ik*  r 
e 


(k,  P) 


With  this  approximation  the  "scattered  wave"  ^(r)  is 
given  by  the  Lippmann  -  Schwinger  equation  (13) 


elk‘ (r  r< )y(rT )0 (£-r* )^(rt ) 

e(kQ*  "?)  -  e  (£,■?)  ±  i€ 


(2.11) 


, 
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The  center  of  mass  momentum  P  and  the  relative 
momentum  It  are  defined  in  terms  of  the  single  particle 
momenta:,  Ic-^  and  Ic^  by 

*=!(*!  - $2) 

and 


The  terms  in  the  energy  denominator  of  equation  (2.1l) 
are  given  explicitly  by 


e  (If,?)  =  E  (|  +  It)  +  E  (f  -  If) 


j-  +  k2  +  U  (|  +  f )  +  U  (|  -  Ic) 


(2.12) 


where  E(k^)  and  U(Ic-j)  represent  the  energy  and  the 
potential  energy  respectively  of  a  particle  with 
momentum  Ic ^ .  We  now  assume  that  the  energy  denominator 
may  be  represented  by  a  quadratic  function  of  k, 
i  .  e . ,  we  put 


e  (fo,5,?)  =  e  -  e  (1?,?) 

A(ko)  +  k2 

_  -  - 

m 


(2.13) 


We  do  not  attempt  to  justify  this  assumption  here.  It 
is  useful  in  that  it  enables  us  to  derive  an  equation 


1 6. 


for  two  particle  scattering  which  is  later  seen  to  be 
identical  with  an  equation  obtained  for  the  nuclear 
many-body  problem  by  the  reference  spectrum 
approximation.  Equation  (2.1l)  may  now  be  written 


ik  •  r 

f (r)  =  e  °  + 


* 


m 


(2t r) 


/dr'  /die 


;ik  r  ^y(r T  )d  (^-r T  )il/(r'  ) 


-  A  (k  )  -  k‘ 
v  o  1 


±  ie 


(2.14) 


The  ±i€  term  in  the  denominator  is  chosen  to  give 
incoming  or  outgoing  scattered  waves.  But  from  equation 
(2.9)  it  follows  that  there  must  be  no  phase  shift  due 
to  the  short  range  potential.  Then  ^(r)  must  be  a 
standing  wave  and  equation  (2.14)  must  be  solved  by 
taking  the  principal  value  of  the  integral.  Then  we 
can  transform  (2.l4)  to  a  differential  equation 

p  ix  ♦  ?  * 

[ v  -  A (kQ ) ]  [^(r)  -  e  ]  =  m  v(r)0(£-r)  ^(r)  .  (2.15) 


The  wave  function  ^( ~r)  may  be  expressed  as  a  sum  of 
partial  waves 


f  (r)  =  7  2 

r  1=0 


(2 i  +  l)il  P£  (cos  9)  u^(r) 


(2.16) 


17. 


The  equation  for  the  radial  wave  function  u^(r)  is 


(2.17) 


where  J.(x)  is  the  spherical  Bessel  function  of 


the  first  kind. 

Since  the  potential  is  infinitely  repulsive 
for  r  <  c  ,  the  wave  function  must  vanish  at  r  =  c. 

From  the  fact  that  there  is  no  phase  shift  we  have  another 
boundary  condition.  These  plus  the  fact  that  the  wave 
function  and  its  first  derivative  must  be  continuous  at 
the  separation  point  are  the  boundary  conditions  on  the 
solution  of  equation  (2.17). 


Ug(°)  =  0 

U^(l’)  =  UgU4  ) 
Ug ’  ( £ ~ )  =  ’  U+ ) 


(2.18) 


(2.19) 


(2.20) 


(2.21) 


For  r  >  £  the  solution  of  equation  (2.17)  which 
satisfies  condition  (2.21)  is  (assuming  A  <  0) 


Uj(r)  -  rj^(kr)  =  B  r  (i/  A  r) 

(1)  ,  , 

where  h^  [x)  is  the  spherical  Hankel  function. 
Hence  in  solving  equation  (2.17)  for  r  <  £  9  the 
boundary  conditions  would  be  equation  (2.18)  and 


dr  {  Vr)  -  rW)3 

"a  ,  H) 

-5—  1  n 
dr  l 

(i ]  A(kQ)  r) } 

Vr)  ■  rh(kor) 

r=| 

r  h^  (j 

J  A(kP  r3 

(2. 

The  solution  must  be  consistent  with  equation  (2.9) 
which  has  the  following  form  in  terms  of  the  partial 


waves , 
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Chapter  3«  The  Nuclear  Many-Body  Problem 


Assuming  that  the  Hamiltonian  contains  just 
the  kinetic  energy  and  a  two-body  interaction,,  we 
have  for  a  system  of  A  nucleons 

A  A 

H  =  2  T  +  i  2  v(r  )  (3.1) 

1=1  1  ^  l,j=l 


Ih  is  the  kinetic  energy  operator  for  a  particle  in 

the  state  i.  In  the  second  term  it  is  to  be  understood 

that  self -interaction  terms  are  excluded  from  the 

summation.  A  local  potential  has  been  assumed 

(r.  .  =  r„  -  r.).  The  Hamiltonian  is  divided  into  two 

parts 

H  =  H  +  H 1 
o 


where 


H  = 
o 


2 

1 


T.  + 


1  2 

2 


v (r.  . )  6 (r.  . 
v  ij  ij 


-  | ) 
ij 


(3-2) 


and 


H  ’  —  — 

2 


2 
i  j 


v  (?ij)  e(Hj  -  rij) 


(3-3) 


The  unperturbed  Hamiltonian  Hq  contains  the  kinetic 
energy  and  the  long  range  attractive  potential.  The 
perturbation  HT  consists  of  the  short  range  potential 


which  is  strongly  repulsive  for  very  small  r 
strongly  attractive.  The  separation  distance 
not  yet  specified. 


and  then 


20. 


.  is 
i  j 


The  total  Hamiltonian  H  defines  the  exact 
wave  vector  |  >  and  the  total  energy  E  of  the 

ground  state  of  the  system 


H|  T  >  =  E  |  T  >  , 
o  o 


(3.4) 


and  H  defines  a  set  of  unperturbed  wave  vectors 
o 

|  0^  >  with  eigenvalues  E^ 


H  |  $  >  =  E  |  0  >  . 

o  n  n  n 


(3-5) 


The  exact  wave  function  |  T  >  must  be  normalized  such 

o 

that 


<  0  |  1  >  =  1  . 
o  o 


(3.6) 


Then  I  T  >  is  expanded  in  terms  of  the  complete 


o 


orthonormal  set  |  $n  >  to  yield 


T  >  =  |  0  >  +  2. 

o  o  n4=0 


0  >  <0  |  H1  |  T  > 

n  n  o 


E  -  E 


n 


(3.7) 


E  =  E  +  <4>  I  H'  I  S'  > 

o  o  o 


and 


(3.8) 
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The  energy  shift,  AE  =  <0  |  HT  |  T  >  ,  is  a  function 

o  o 

of  the  separation  distance  £.  .  and  £.  .  can  vary  from 

zero  to  infinity  depending  on  the  momenta  k^  and  kj 

of  the  two  interacting  particles  as  shall  be  shown 

later.  At  this  point  it  is  noted  that  if  £.  .  =  c, 

i  J 

where  c  is  the  radius  of  the  region  of  strong 

repulsion  (the  hard  core  radius),  then  AE  (£  =  c)  >0. 

On  the  other  hand  if  £  is  very  large,  so  that  the  long 

range  potential  is  very  weak,  then  the  major  contribution 

to  E  will  be  the  kinetic  energy  and  E  will  be 

positive.  Since  we  demand  that  the  system  be  bound, 

i.e.  E  <0,  this  means  that  for  sufficiently  large  £  , 

AE  must  be  negative.  Thus  if  the  energy  shift  varies 

continuously  with  it  is  possible  to  make  a 

non-trivial  choice  of  £.  .  in  such  a  way  that  the  energy 

J 

shift  vanishes. 


For  each  pair  of  particles  a  wave  vector 

T.  ,>  is  now  defined  such  that  |  T.  .  >  contains 
i  j  iJ 

correlations  of  all  pairs  except  ij 


T.  .>  =  |  $  >  +  z 

ij  o 


|  <D  >  <0  |  H*  -  v.  .  |  T  > 

1  n  n  l ,]  o 


n+0 


E  -  E 


n 


(3-9) 


where  v.  .  =  v(r.  . )  e(£L  ,  -  r.  .)  is  the  short  range 
1J  1J  IJ  IJ 

s 

part  of  the  potential.  The  operator  G.  .  is  defined  by 

“*■  J 
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s 

G.  .  .> 

ij  ij 


s 

V.  . 

IJ 


*o> 


(3.10) 


Equations  (3. 10)  and  then  (3*9)  may  be  substituted  in 
(3-7)  to  yield 


Y  > 
o 


0  >  + 
o 


2 

n^=0 


$  >  <0 
n  n 


1 

2 


l 


GS 


il 


|  $  > 


o 


E 


-  E 


n 


+ 


(3.11) 

where  the  remaining  terms  are  contributions  from 

three,  four,  and  many  particle  excitations. 

Neglecting  these  higher  order  excitations  is  equivalent 

to  assuming  that  |  ¥.  .>e£  |  $  >  .  Equations  (3*7) 

ij  o 

and  (3*9)  give 


Y  >  =  I  1.  .>  + 
o  ij  n*=0 


0  >  <0  |  v.  .  |  W  > 

n  n  '  ij  o 


(3.12) 


E  -  E 


n 


We  let  v^j  operate  on  this  equation  and  substitute 

for  v.  .  I  ¥  >  from  equation  (3.10) 
ij  o 


G?  .  |  I.  .  > 
ij  ij 


s 

v.  . 

ij 


v?  .  U  >  <*  I  G?  .  I  .> 
¥  >  +  z  1J  n  n  1 J  iJ 
i  j  n*0 


E  -  E 


n 


(3.13) 
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Then  with  the  approximation  I  T.  .>  |  $  >.,  one 

ij  -  o 

g 

obtains  the  integral  equation  for  Gj_j 


Q  <$  |  v?  .  |  0  >  <$  |  G?  .  I  $>  > 

i  s  |  .  x  .  v  m  ij  1  n  n'ij'o 

<3>  G.  .  0  >  =  <$  v.  .  3>  >  +  ^ d - 

m  ij  1  o  m  ij  1  o  n*0 


E  -  E 


n 


The  separation  distance  £.  .  is  now  specified 

J 

by  imposing  the  condition 


(3.14) 


<  0  |  G7  .  |  $  >  =  0 

o'  ij  o 


(3-15) 


Then  with  the  approximation  |  ^  .  >  =  |  ,  the 

energy  shift  is 


<*o  I  H'  I  V  =  I  U  <$o  1  vlj  1  V 


I  ij  <*o  I  GiJ  I  *ij> 


* 1  <h  i  G?j  i  v 


=  0 


and  the  total  energy  of  the  system  is 
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E  =  E 


o 


=  2 
1 


<0 


o 


T 


2 


v 


i 

ij 


V 


2 

i 


E. 

l 


|  .2.  <$ 


o 


$  > 
o 


(3.16) 


where  v.  .  =  v(r  )  0(r .  .  -  £.  .)  is  the  long  range 

J  J  J  1  J 

part  of  the  potential  and  E^  is  the  single  particle 
energy 


E. 

1 


<$ 


o 


T,  |  0  >  +  .2.  <0  |  v.  .  |0  > 
1  O  J*1  o'  1J  1  o 


(3-1?) 


It  is  now  assumed  that  equation  (3.14)  can  be 
reduced  to  a  two-particle  scattering  equation.  Consider 
I  <$n>  to  be  specified  by  the  single  particle  wave 
functions  I  (i.e.  a  Slater  determinant).  Then 

g 

the  matrix  elements  of  G.  .  depend  only  on  the  single 

J 

particle  wave  functions  of  particles  i  and  j  and 
the  diagonal  elements  vanish  from  equation  (3.15). 

Thus  in  (3.14)  one  need  only  sum  over  all  possible 
excitations  of  particles  i  and  j.  The  energy 
denominator  will  be  the  difference  between  the  energies 
of  particles  i  and  j  in  the  "ground  state"  |  ij> 
and  the  excited  state  |  mn> 


25. 


T  T 

<i  J 


r  r 


Gij  I  I  v^j  I  i  j>  + 


mn+i  j 


!  T 

<i  J 


s 

V.  . 

IJ 


mn>  <mn 


G? 


IJ 


I  ij> 


ij 


-  e 


mn 


(3.18) 


The  separation  distance  £.  .  must  now  be  such  that 

i  J 


<  ij  |  G  |  ij>  =  0 


(3.19) 


is  satisfied.  If  the  separation  distance  is  known 
then  the  unperturbed  Hamiltonian  Hq  is  exactly 
determined.  Since  Hq  contains  only  the  long  range 
potential  which  is  non-singular  one  could  use  the 


Hart ree 


Fock  method 


(i4) 


to  calculate  the  single 


particle  wave  function  |  cp.>  and  thus  determine  the 

energy  of  the  system  which  will  be  the  eigenvalue  of 

H  .  This  is  impractical  however  since  a  very  accurate 

knowledge  of  is  required  and  to  solve  equations 

(3.18)  and  (3.19)  for  4-  •  one  must  presuppose  a 

J 


knowledge  of  the  wave  function  |  qu>  .  One  would  have 
to  obtain  a  self-consistent  solution  through  an  iteration 
and  the  iteration  could  not  be  expected  to  converge 
very  rapidly  because  the  Hartree  -  Fock  solution  for  |  qu> 


depends  so  critically  on  £. 


IJ 


An  alternate  method  of 


attacking  the  problem  would  be  to  use  a  perturbation 


expansion  of  H  .  But  if  there  is  a  strong  tensor  force 
the  expansion  may  not  converge  rapidly  enough and 
the  calculation  would  have  to  be  carried  out  to  third 
order  at  least. 

At  this  point  the  plane-wave  approximation 
for  |  cp.  >  is  introduced,  i.e.,  it  is  assumed  that 

ilc.  •  r . 

|  cp .  >  is  a  momentum  state  I  Ic .  >  and  that  <r.  I  cp.>  =  e 
What  this  means  is  that  the  ideal  fermi  gas  is  chosen 
as  the  point  of  departure  and  the  wave  function  |  0q> 
is  approximated  by  the  ground  state  wave  function  of 
a  fermi  gas,  i.e.,  a  determinental  wave  function  in 
which  all  plane  wave  states  with  wave  number  less  than 
or  equal  to  the  fermi  momentum  k^  are  occupied  and 
all  other  states  are  unoccupied.  The  validity  of  this 
approximation  has  been  demonstrated  by  Bethe  (l 6)  for 
a  well-behaved  central  force  and  by  Moszkowski  and  Scott  (9) 
for  a  hard  core  plus  an  attractive  central  potential. 

Thus  equations  (3*18)  and  (3-19)  are  written 


<Jt.  ,Ic .  |  G?  .  |  If.  Ic  .> 

i  J  ij  i  j 


<£.  I  VS  .  I  it,  Jc  .> 

1  J  '  IJ  i  J 


f  t 


+ 


2 

mn*i  j 


<gi 


v, 


•  • 

2JL 


Ic  1?  >  <1?  Tc  |  Gf  .  |  Ic.  lt.> 
m  n  m  n  i  ,i  1  l  j 


e  (k±,kj)  e  (kmI<n) 


(3-20) 
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and 


(3.21) 


The  prime  in  equation  (3.20)  indicates  that  the  summation 
is  only  to  be  taken  over  those  excited  states  which  are 
permitted  by  the  exclusion  principle.  An  operator  Q 
is  inserted  to  take  care  of  the  exclusion  principle 
where  Q  =  1  if  k  and  k  are  greater  than  k^  and 
Q  =  0  otherwise.  It  is  also  convenient  at  this  time 
to  introduce  the  relative  and  center  of  mass  momenta 
by  putting 


(3.22) 


and 


(3-23) 


"?  =  k .  +  Ic .  . 

1  J 


Then  since  the  total  momentum  of  the  two  interacting 


particles  must  be  conserved,  will  merely  be  a 
parameter  of  the  integral  equation  for  Gs .  Hence  the 
equations  (3.20)  and  (3-2l)  are  written 


<lx  |  gs  |  =  <17  |  vs  |  o>  +  2 


<£  |  Vs  |  Ic T  >Q(lq  ?,  kp),<£T  |  Gs  |kQ,  P> 


and 


e  (£q,  £  ,  T) 

(3-24) 


(3.25) 
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The  energy  denominator  is 


e  (I?0,  If',  T)  =  E  (It)  +  E(It.)  -  E(I?m)  -  E(Icn) 

=  E(|+  I?0)  +  E(|  -  lfo)  -  E(|  +  If')  -  E ( 


(3-26) 


where 


represents  the  single  particle  energy.  For 


particles  in  the  fermi  sea  equation  (3. 17)  becomes 

2 

l 


k. 

e  (ki )  = 


+  Z  <k  , k 
1  J 

J  1 


v: .  k. j k . -  k  . , k. >  .  (3.27) 
ij'i  j  j  1  *  '  ’ 


The  single  particle  energy  for  excited  states  will 
be  defined  later.  It  is  also  possible  to  obtain  an 
integral  equation  for  the  perturbed  two  particle  wave 
function  by  putting  |  I  i n  equation  (3.12) 

and  assuming  E  -  E^  depends  only  on  the  excitation 
energies  of  the  two  particles.  The  result  is 


f(’£o,7)  > 


I  )<f'  |  vs  |  f( k  ,?)> 

$,?>+!.  - - - 2 - 

0  If’  e  (Jco,  V ’ ,  T) 

(3.28) 


g 

The  above  equations  can  be  solved  for  i,  G  ,  and  f . 

The  equations  must  be  solved  self -consistently  for  an 
exact  solution.  The  method  of  obtaining  the  solution  will 


t  p^Ioj 
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be  discussed  in  the  next  chapter.  It  should  be  noted 
here  that  4  will  be  a  function  of  the  relative  momentum 
kQ  and  will  also  depend  on  the  parameters  P  and  k  . 

g 

The  purpose  of  introducing  the  matrix  G 
and  defining  the  separation  distance  £,  is  to 
facilitate  calculation  of  the  Brueckner  reaction  matrix 
G  which  is  defined  by  an  integral  equation  similar  to 
equation  (3.24),  except  that  the  complete  potential  is 
used.  In  operator  form  the  equation  for  G  is 

G  =  v  +  v  —  G  (3-29) 

e 

The  matrix  elements  of  G  are  also  given  by 

<k  |  G  |  £q,P>  =  <Tc  |  v  |  ^N( (3.30) 

where  |  ^N(Tcq^)>  is  the  wave  function  for  two  nucleons 
interacting  through  the  potential  v(r).  The  integral 
equation  for  |  is 

Tc*>  Q<£  I  v  |  ^N(l Tq,1?)> 

e(£Q,£T, 

(3.31) 

For  brevity,  the  parameter  ?  will  be  omitted 

From  the  definition  of  Gs,  it  is  seen  that 


hereafter. 


30. 


equation  (3.25)  may  be  written 


<  |  vs  |  ^(Ic  )>  =  0  . 

o  1  v  o 

Thus  the  diagonal  elements  of  G  (equation  (3.30)  )  are 
given  by 

<  T?o  |  G  |  Tfo>  =  <£o  |  v*  |  Ai?o)>  ,  (3.32) 


with  the  approximation  \f>  =  |  being  made.  A 

comparison  of  equations  (3.28)  and  (3 • 31 )  shows  that 

o 

this  is  a  good  approximation  if  v  is  weak.  Now  it  is 
well  known that  at  large  distances  (i.e.  r  >  £  ) 
the  free  particle  wave  function  is  a  good  approximation 
to  the  two  particle  wave  function  in  nuclear  matter. 

Thus  to  first  order  we  have 

<  T?o  |  G  I  Tf0>  =  <J?0  I  wl  |  1?0>  (3-33) 


More  exactly,  a  comparison  of  the  integral  equations  for 
G  and  Gs  shows  that 


&  , 
+  v  + 


(Q  as 

v  e 


+ 


l  ,  J  Q  ns  J  Q  i 

v  +  v  —  G  +  v  —  V  + 
e  e 


.  •  (3.34) 


where  the  remaining  terms  are  of  third  and  higher 
order  in  Gs  and  v^.  Thus  to  first  order  G  =  Gs  +  v^ 
and  (3-33)  gives  the  diagonal  elements  of  G.  The  importance 
of  this  result  is  that  to  first  order  the  diagonal  elements 
of  the  reaction  matrix  can  be  calculated  from  the  long-range 


31. 


potential  alone.  Of  course  the  short  range  potential 
must  be  known  for  an  exact  calculation  of  the  separation 
distance  £  but,  as  will  be  shown  in  chapter  4,  a  first 
order  calculation  of  the  separation  distance  can  be  based 
on  the  known  two-nucleon  scattering  phase  shifts  and 
the  outer  potential  only.  Thus  it  is  possible  to  make 
an  approximate  calculation  of  the  reaction  matrix 
without  knowing  the  short  range  potential. 

We  now  wish  to  use  equation  (3-33)  to  express 
the  reaction  matrix  as  an  integral  of  the  long  range 
potential.  To  obtain  the  correct  form,  one  must  at  this 
point  consider  the  spin  and  isospin  states  of  the  two 
nucleons,  which  are  required  for  the  complete  specification 
of  the  dynamical  state  of  the  two-particle  system.  In 
order  to  have  the  reaction  matrix  G  depend  only  on 
the  momentum  states  of  the  interacting  pair,  G  is  now 
specified  as  an  average  over  the  spins  and  isospins  of 
the  interacting  pair.  Writing  the  exchange  term 
explicitly  we  have 

<  |  G  |  Tc  >  =  [  <  ij  |  G  |  ij  -  ji  >  )  (3.35) 

0  °  spin, isospin  average 

Now  if  the  spin  and  isospin  are  ignored  and  the 
exchange  term  is  omitted,  one  can  use  the  partial  wave 
expansion  of  the  plane  wave. 
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i  Tr  •  7?  00  n 

e  =2  (2£  +  l)  i  (cos  0)  (kr)  , 

H/— 0 

to  express  <  k  |  G  |  k>  in  the  form  (for  zero  center  of 
mass  motion) 

00 

4tt  Z  (2£  +  l)  f  (kr)  v(r)  r2  dr  , 

£  5(k) 

where  the  normalization  factor  (inverse  of  the  volume) 

has  been  omitted  for  convenience  in  passing  to  the 

limit  of  infinite  volume.  The  correct  form  to  use  will 

be  the  above  with  statistical  weights  c^  included 

to  correct  for  the  exchange  term  and  the  spin,,  isospin 

averaging.  The  statistical  weights  are  deduced  as 

follows:  for  given  momentum  states.  Ic.  and  Ic..  the 

i  J 

spin,  isospin  states  of  the  two  particles  will  differ 
three  fourths  of  the  time  allowing  all  values  of  Z  to 
contribute  to  the  direct  term  and  causing  the  contribution 
of  the  exchange  term  to  vanish  since  the  spin  states  are 
orthogonal.  One  fourth  of  the  time  the  spin,  isospin 
states  of  the  two  particles  coincide  and  the  exchange 
term  cancels  the  direct  term  for  even  i  contributions 
and  gives  a  factor  two  for  odd  &  contributions.  Thus 
with  the  plane-wave  approximation,  the  diagonal  elements 
of  G  as  defined  by  equation  (3.35)  are 


33. 


II  00  00  Q  Q 

<k  |G|k>  =  4tt  2  c.  (2,0  +  l)  f  j  „  (kr)v(r)r^dr 

i=0  ^  |(k) 


(3-36) 


3  3 

where  ^  for  even  X  and  ^  for  odd  X. 


This  expression  is  valid  for  a  central  potential  which 
is  independent  of  the  spin  and  isospin  of  the  interacting 
particles.  A  more  realistic  inter-nucleon  potential 
depends  on  the  total  angular  momentum,  orbital  angular 
momentum,  spin,  and  isospin  of  the  interacting  nucleons. 
Denote  such  a  potential  by  v ^  j  (r)  where  J,  Z,  S,  and  T 
denote  the  above-mentioned  parameters  in  that  order 
according  to  the  usual  convention,  i.e.  J(J  +  l)  is 
the  eigenvalue  of  J  ,  etc.  For  such  a  potential  the 
reaction  matrix  is  given  by 


<Ck  |  G  |  =  8tt  y7~  2 

16  odd  Z 


(2 Z  +  l)  /  j ^(kr)  v 

£ 


+  T7'  2 

D  even  Z 


00 

(2.0+1)  /  (kr)  (r)  r2dr 

£  3 


i  2 


2 

even  Z 


odd  ^ 


(3.36) 
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This  equation  follows  from  the  observation  that  the 
exchange  term  cancels  the  direct  term  for  states  of 
even  parity  (£  +  S  +  T  even)  and  gives  a  factor  two 
for  states  of  odd  parity  (£  +  S  +  T  odd).  Thus  only 
states  of  odd  parity  are  included  and  a  factor  two  is 
inserted  outside  the  square  bracket.  The  T  =  1  terms 
are  multiplied  by  three  to  account  for  the  isospin 
degeneracy  and  the  spin  degeneracy  is  taken  care  of  by 
the  summation  over  J  in  the  S  =  1  terms.  Finally  the 
spin,  isospin  average  is  obtained  by  dividing  throughout 
by  sixteen  which  is  the  number  of  possible  spin-isospin 
states  of  two  nucleons.  Equation  (3-36)  can  be 
simplified  by  defining  two  potentials  v^(r)  and 
v^(r)  by  the  following  relations 

+  1)  vt (r)  =  3(2i  +  1)  v°'kr)  +  2  (2J  +  l)  v*’°(r) 

'  3  J  3 

[2H  +  1)  v'(r)  =  (?-i  +  1)  v°’°(r)  +  3  j  (2J  +  l)  v^(r) 
Then  one  obtains 

00  o 

<k  |  G  |  k>  =  |  [  2  (2i  +  1)  /  (kr)  vj(r)  r  dr 

odd  l  ‘ 

00  „ 

+  2  (2j0  +  1)  /  (kr)  v*(r)  r  dr] 

even  ^  £ 


(3-37) 

(3.38) 

(3-39) 
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+ 

If  the  potentials  v ^  and  are  independent  of 
this  can  be  further  simplified  to 


<k  |  G  |  k>  =  £  /  [ ( v+ ( r )  +  v  (r))  +  (v+(r) 

hk) 


v  (r))j  (2kr)]  kdr 
(3-40) 


We  now  consider  the  energy  denominator  which 

g 

occurs  in  the  integral  equations  for  G  and  G  .  From 
equation  (3.2 6)  one  can  write 


e(CoJlf‘,?)  =  k^  -  k'2  +  U(|+  i?0)  +  U(£  -  -  U(|+  If')  - 

- ► 

-  U(|  -  Jc')  (3.41) 

where  Ufk^)  is  the  single  particle  potential.  This 
function  is  always  defined  in  the  literature  by 


U(k. )  =  2  <  k  |  G  |  k  > 

kj-<kp 


(3-42) 


wher^e  k  =  i  (k.  -  k.).  This  can  be  written  differently 

2  v  l  j 

to  show  the  exchange  term  explicitly,  thus 


U(Tc.  )  =  2  <£.  |  G  |  -  Icj ,\z±>  (3*^3) 


^  k  <^k  ^  ^ 

j  N  F 
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The  meaning  of  this  function  becomes  more  apparent 
when  we  take  the  matrix  element  of  equation  (3.1) 
between  the  wave  vector  of  the  system  I  TN>  and  the 
unperturbed  eigenvector  |  0Q>  ,,  keeping  in  mind  that 
<0o  I  =  1.  One  obtains 


E  = 


Z 

i 


<0 


o 


To 

l 


iN> 


+ 


1_ 

2 


2 

•  • 

li  J 


<0 


o 


V, 


ij 


tn> 


(3.44) 


Now  I  0  >  has  been  assumed  to  be  a  Slater  determinant 
o 

of  plane  waves  and  |  is  the  anti symmetrized  nuclear 

many-body  wave  vector  with  two  body  correlations  given 

by  equation  ( 3 • 31 ) *  Thus  I  $Q>  is  an  eigenvector  of 

1  2 

T.  with  eigenvalue  77k.  and  for  the  potential  energy 

1  C-  1 

contributions  from  v . .  one  can  integrate  over  all 

—  J 

co-ordinates  except  r.  and  r.  .  In  this  way 

--  j 

equation  (3.44)  is  reduced  to 


k. 


E  = 


Z 

l 


2 


4  2  <E  ,£.  |  V  |  fN(i?,n  >  (3.45) 

i  j  j  J 


where  k  and  ?  are  the  relative  and  center  of  mass 
momenta  of  particles  i  and  j  as  defined  by  equations 
(3.22)  and  (3.23).  Then  with  equation  (3-30)  this 
expression  can  be  reduced  to 
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E 


1 

2 


2  <  |G  |  17  > 

i  j 


(3.46) 


=  2  [E(ic.)  -  I  u(£.)] 

i 


(3.^7) 


This  definition  (equation  (3.42))  of  the  single  particle 
energies  is  more  exact  than  the  expression  used 
previously  (equation  (3.17))  and  from  equation  (3.34) 
it  follows  that  to  first  order  the  definitions  are 
equivalent.  This  energy  spectrum  is  not  intended  to 
give  the  actual  energy  required  to  remove  a  nucleon  of 
momentum  k.  ^  k^  from  the  nucleus,  rather  it  is 
designed  so  as  to  yield  an  accurate  value  for  the  total 
energy  of  the  nucleus  from  equation (3 • 47) .  For  the 
energies  E(it  )  of  virtual  states  above  the  Fermi 
sea  we  put 


E(£  ) 

v  nr 


k 


m 


+ 


k  <k 
nx  F 


<  k  ,  It 
m  n 


G 


m  n 


it  *  It  > 
n  m 


(3 


Once  again  this  has  nothing  to  do  with  the  actual  energy 
of  a  real  nucleon  in  an  excited  state.  The  definitions 
of  E(k.  )  and  E(km)  are  adopted  for  internal  consistency 
in  the  calculations.  The  energy  denominator  is  the 
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difference  in  the  total  energy  of  the  nucleus  in 
the  ground  state  and  in  the  excited  or  virtual  state  and 
the  single  particle  energies  are  so  defined  that  the 
energy  denominator  will  be  given  by  equation  (3.4l). 


To  find  the  average  potential  of  a  particle 

with  momentum  k^,  therefore,  one  calculates  the 

average  value  of  G  when  k^  is  fixed  and  kj  takes 

on  all  values  in  the  Fermi  sea.  The  summation  over  the 

states  |  !?.>  reduces  to  a  multiplication  of  this 
J 

result  by  the  density  p  (instead  of  the  number  of 
states,  since  the  reaction  matrix  elements  have 
already  been  multiplied  by  the  volume  for  convenience 
is  passing  to  the  limit  of  infinite  volume)  where 
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4t t  k 
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2  k 
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(2t r) 
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3tt 
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(3.49) 


A  factor  of  4  is 
isospin  degeneracy, 
thus  p 

u(ki)  =  - 


included  because  of  the  spin  - 
The  single  particle  potential 
f  <k  |  G  |  k>  d  ic. 

J&F  1 

/  d  ir 


is 


(3.50) 
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On  substituting  equation  (3.40)  for  G,  this  becomes 


U(k±) 


1 


/ 


8"  kJ<kF 


dk . 
J 


00 


I 


[  V+  +  v' 


(k) 


+  (v+-v  )  j  (2  kr)]  r^dr 


(3-51) 


where  k  =  ^  I  It  —  Tc.  I  . 

2  1  i  J 


It  is  now  possible  to  discuss  the  average 
binding  energy  per  nucleon.  For  this  we  need  the 
average  potential  energy  of  the  nucleons  in  the  ground 
state  which  can  be  written  in  terms  of  U(k^) 


U 


/ 


k±^k 


U(k.  )  d  It 

v  l '  l 


/ 

k^kp 


d  k. 
i 


(3.52) 


or  directly  as  an  integral  of  G 

kp 

U  =  p f  dk  P(k,  kp)  <  k  |  G  |  k> 
0 


(3.53) 


where 


P(k, 


kp) 


24  k‘ 

1  r  3 


(  1  - 
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2  k 


+ 


k 
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2  kn 
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(3.54) 


is  the  normalized  probability  distribution  of  relative 
momentum  for  two  particles  in  the  Fermi  sea.  The 
binding  energy  per  nucleon  is 
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W (kp)  =  T  +  i  U 

=  fo  kF2  +  |  U  (3.55) 

For  low  densities.,  instead  of  equation  (3*53)  for  U, 

one  can  use  the  reaction  matrix  evaluated  at  the 

average  relative  momentum  k  multiplied  by  the  density,  i.e., 

k  3  00 

U  =  /  [  v+  +  v"  +  (v+-v”)  jQ(2  kr)]  r2dr  .  (3-56) 

£(k, kp) 

The  average  relative  momentum  is  a  function  of  k 

F 

which  depends  on  the  form  of  <k  |  G  |  k>  .  Usually 
one  uses  k=  v/o  .3  kp  which  is  exact  if  <k  |  G  |  k> 
is  a  simple  quadratic  in  k.  With  the  replacement 


I  (k  (kp) ,  kF)  =  %  (kj 


(3.57) 


the  binding  energy  is 

o 


w  (k„)  =  ^  k y  + 


k. 


00 


/  [v+  +  v 

1  (kj 


10  F  12t r 


+  (v+-v  )jQ(2  kr)  ]r  dr 


(3.58) 


It  is  now  possible  to  discuss  the  problem  of 

the  existence  of  the  minimum  of  W  (kp). 


saturation,  i.e.. 
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in  terms  of  conditions  on  the  derivatives  of  f(kp)  . 
To  simplify  the  equations,  it  is  assumed  that  the 
term  (v+  -  v~)jQ(2  kr)  does  not  contribute  to  the 
reaction  matrix.  Then,  putting 


v+  (r)  +  v  (r)  =  V(r) 


we  have  ^ 

k  d  00 

W(kp)  =  kF2  +  if?  J  V(r)  r  dr  (3’ 

?  (kp) 

The  saturation  condition  states  that  at  the  equilibrium 
density,  i . e .  at  kp  =  kpN,  W(kp)  should  be  a  minimum. 
Thus  at  the  nuclear  fermi  momentum,  kpN  ,  W(kp) 
should  satisfy 


d  W 

d  kp 


0 


>  0 

A  straightforward  substitution  of  equation  (3-59)  shows 
that  in  terms  of  the  separation  distance  £,  (kp)  these 


d2  W 

d  kp 
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conditions  are 
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(3.61) 


These  equations  do  not  help  in  the  determination 
of  £  (k^)  but  they  would  be  very  useful  as  a  means 
of  checking  the  theory.  It  would  be  particularly 
interesting  to  evaluate  the  separation  distance  to 
first  order  using  only  the  long  range  part  of  the 
potential  and  see  if  reasonable  predictions  can  be  made 
without  any  knowledge  of  the  form  of  the  short  range 
potential . 

It  is  also  of  interest  to  compare  the  binding 

energy  with  the  energy  of  a  particle  at  the  top  of  the 

Fermi  sea.  According  to  a  theorem  of  Hugenholtz  and 
( 1  7 ) 

Van  Hove v  '  these  should  be  equal.  The  total  energy  of 


the  system  is  given  by  equation  (3-46)  which  may  be 
written  as  an  integral  in  the  form 


4 


3 


E  =  /  [T(k  )  +  iu(k  )  ]  d  k  (3.62) 

^  k.«kp  121 
l  F 


where  ^  is  the  volume  of  the  system.  Then  the 
energy  per  unit  volume  is 

kF 

E =  -|  /  [k  +  u  (k  )  ]  k  2  d  k  .  (3.63) 

7 r  0 


The  binding  energy  may  be  obtained  by  differentiating 
this  with  respect  to  the  density  of  nucleons^1^,  i.e. 

d  E  k-p 

W  =  - -  =  — /  [k  2  +  U  (k  )  ]  k  2  d  k  (3-64) 

dp  2  kp  a  kF  0  1  1  1  1 

Then  if  equation  (3.5l)  is  substituted  for  U  (k^)  this 
becomes 

W  =  T  (kp)  +  U  (kp) 


64 


I  f  i2  [v++v  +(v+-v")  J0(2k4)  jd^dltj 

k^ } kj^kp 


(3.65) 
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where  the  potentials  are  evaluated  at  r  =  £(k,  kp) 
and  k  =  75-  |  Tc ^  -  ic j  |  .  The  difference  between  W  and 
E(kp)  is  called  the  rearrangement  energy.  It  should 
be  noted  that  the  first  order  calculation  of  £,  from 
the  long  range  potential  only  gives  a  result  which  depends 
only  on  relative  momentum  k  and  for  this  i 3  there 
will  be  no  rearrangement  energy.  If  £(k*  k^)  is 
accurately  evaluated,,  there  may  be  a  non-vanishing 
rearrangement  energy  which  can  be  attributed  to  the 
approximation  used  in  the  expression  for  U(k^)  ,  i.e., 
the  neglect  of  the  higher  order  terms  in  the  expansion 
of  the  reaction  matrix  (equation  (3*34)  ).  The  value  of 
the  rearrangement  energy  might  then  be  expected  to 
provide  some  idea  of  the  relative  importance  of  these 
higher  order  terms. 
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Chapter  4.  The  Calculation  of  the  Separation  Distance 


It  has  been  shown  in  the  previous  chapter 
that  to  a  good  approximation  the  diagonal  elements  of 
the  reaction  matrix  are  given  by  the  first  Born 
approximation  of  the  long  range  potential.  With  the 
problem  formulated  in  this  way,  the  solution  of  the 
Brueckner  equations  is  reduced  to  the  problem  of 
calculating  the  separation  distance  £  in  a  self- 
consistent  way. 


A  first  order  calculation  can  be  based  on  the 
knowledge  of  the  two-nucleon  scattering  phase  shifts 
and  the  long  range  part  of  the  potential.  The  first 
order  approximation,  £  ,  is  obtained  by  letting  the 
interaction  energy  in  the  energy  denominator  tend  to 
zero.  Also  the  center  of  mass  of  the  two  interacting 
particles  is  assumed  to  be  at  rest  and  the  exclusion 

operator,  Q,  is  replaced  by  the  principal  value  operator. 

Bs  s 

Then  G  ,  the  first  order  approximation  to  G  satisfies 

the  equation 


(4.1) 
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F  o  o 

where  e  (k^k1)  =  kQ  -  k*  and  P  stands  for 
Cauchy Ts  principal  value.  This  is  just  the  equation 
for  the  reaction  matrix  of  scattering  theory  with  the 

g 

potential  v  and  the  solution  is  given  in  terms  of 
the  scattering  phase  shifts  due  to  the  short  range 
potential  by 


<k  |  GFs  |  k>  = 


—  47 T 
k 


f  +  1]  ci  5i  (k’  U 


(4.2) 


(From  the  theory  of  scattering  it  is  expected  that 
Fs 

G  should  be  the  reaction  matrix  of  scattering  theory 

g 

(for  potential  v  )  which  is  given  in  terms  of  tan  5 


rather  that 


r 


However  Fukuda  and  Newton 


(18) 


H 


have 


shown  that  the  limit  of  G  ,  as  the  interaction  term 
in  the  energy  denominator  tends  to  zero,  depends  on 
as  shown  above.) 


For  a  spin  and  isospin-dependent  potential,  v 


S’T(r), 


we  have 


<k  I  G 


Fs 


0,1 


k>  = 


-  h  {  z  42i  +  x)  [3  6*  ^ +  36i  (k’i>)] 

even  £ 


1,0 


+  Z  (2£  +  e)  [8°'°(k,|  )  +  (k,e  )]) 

odd  l  HOI 

(4.3) 
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The  phase  function  b^jF(k,r)  is  the  phase  shift  due 
to  the  short  range  potential  (cutoff  at  r)  and 
satisfies  the  differential  equation  (19) 


d  6i'T(k*r) 


dr 


=  -kr2vS*T(r)  [cos  5^T(k,  r)  j ^  (kr)  - 


-  sin  b^T(k,r)  n^(kr)] 


(4.4) 


with  the  boundary  condition 
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S,T 

i 


00  >  =  6®,T(k) 


(4.5) 


S  T 

where  6^*  (k)  are  the  known  scattering  phase  shifts 
produced  by  the  complete  potential  v  ’  (r).  Now  the 
first  order  separation  distance,,  £  (k),  will  be 
specified  by  the  condition  that  the  diagonal  elements 
of  GFs  vanish. 


<  k  |  GFs  I  k>  =  0  . 


(4.6) 


The  solution  is  found  by  integrating  the  first  few 
partial  phase  equations  (e.g.  &  =  0,  1,  and  2)  inwards 
from  infinity  to  a  point  £ where 


' 
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2  {21  +  1)  [3  d£  (k,  io)  +  3  by  *  (k,  io)  ] 

even  H 

+2  {21  +  1)  [6^  (k,  ?Q)  +  9  b’  (k,  iQ)]  =  0 

odd  & 

(4.7) 

The  solution  ?0(k)  can  be  used  to  compute  approximate 
values  for  the  diagonal  elements  of  the  reaction  matrix 
and  the  single  particle  potential  energy  from  equations 
(3*40)  and  (3. 51 ) •  Thus  first  order  results  can  be 
obtained  which  do  not  depend  on  the  knowledge  of  the 
short  range  potential. 

The  first  order  results  for  the  single 
particle  energies  can  be  used  as  a  starting  point  for  a 
self-consistent  calculation  of  i.  The  energy  denominator 
is  then  given  by  equation  (3*^1 )  and  a  more  accurate 
value  for  the  separation  distance,  i,  can,  in  principle, 
be  obtained  by  solving  equation  (3.24)  subject  to  the 
condition  <  k  |  Gs  |  k  >  =  0.  The  calculations  are 
easier  if  we  use  equation  (3.28).  For  zero  center  of 
mass  motion  the  wave  function,  in  co-ordinate  space 
satisfies  the  integral  equation 
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=  elk*r  +  — f  d  ?»  /  £ 


lit*  •  (r-r 1  )  s 


(27 r) 
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v  (rT )#(k, r 1 ) 


kl>k. 


e(k,  kl) 


(4.8) 


From  equation  (3.25),  it  is  seen  that  the  solution  must 
satisfy  the  condition 

/  dr  e"ik  r  vs(r)  r)  =  0  (4.9) 


Then  from  the  first  order  separation 
distance  i*o,  a  first  order  expression  for  G,  U(k1)J 
and  e  can  be  obtained  from  equations  (3*40),  ( 3 • 51 ) 
and  (3.4l).  With  the  new  energy  denominator  a  better 
approximation  for  the  separation  distance  can 

be  found  from  equations  (4.8)  and  (4.9)  and  this  leads 
to  an  even  more  accurate  expression  for  the  energy 
denominator.  The  iteration  is  continued  until  a 
self  consistent  solution  is  obtained.  In  practice  this 
is  a  complicated  method  and  its  only  advantage  over  a 
direct  solution  of  the  Brueckner  equations  is  that  the 
iteration  is  expected  to  converge  quite  rapidly  because 
the  first  order  results  are  close  to  the  final  results. 
Also  it  is  difficult  to  isolate  the  contributions  of 
different  physical  aspects  of  the  problem  such  as  the 
effects  of  the  exclusion  principle  and  the  center  of 


mass  motion. 
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The  reference  spectrum  method  of  Bethe  et 
is  now  introduced  to  simplify  the  mathematics  and  to 
determine  the  contributions  of  the  above-mentioned 
physical  effects.  Equations  (4.8)  and  (4.9)  are  solved 
under  the  following  conditions  which  greatly  simplify 
the  mathematics: 

I.  the  region  of  integration  over  Ec*  is  extended  to 
include  the  Fermi  sphere  and  the  principal  value 
of  the  integral  is  used; 

II.  it  is  assumed  that  P  =  0,  i.e.,  the  center  of 

mass  of  the  two  interacting  particles  is  at  rest; 

III.  the  effective  mass  approximation  is  used  for  the 
energy  denominator,  i.e.,  e(k,  kT)  is  fit  with 
a  quadratic  form  as  a  function  of  k* 

eR(k,  k,>=-720)/  k'2  (4.10) 

m 

Under  these  conditions  equation  (4.8)  can  be  reduced  to 
a  second  order  differential  equation.  The  possible 
singularity  of  the  integral  equation,  which  is  discussed 
in  the  next  chapter,  is  not  considered  here.  The 
resulting  differential  equation  is 
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[v2  -  72(k)]  [*R(t?)  -  elk' r]  =  m*  vs(?)  fR(r)  (4.11) 

In  the  case  of  a  central  potential  we  can  use  the 

"D  — > 

partial  wave  decomposition  of  f  {r) , 

00 

^R(r)  =  ^  2  (2i  +  l)  i^  P y  (cos  0)  (r)  ,  (4.12) 

i=0 

to  reduce  equation  (4.1l)  to  a  set  of  total  differential 
equations 

2 

[  1  ^  -  72]  [Ug(r)  "  r  J^(kr)]  =  m*vs(r)u,(r)  (4.13) 

Equations  (4.1l)  and  (4.13)  are  just  the  two-particle 
scattering  equations  obtained  previously  (equations 
(2.15)  and  (2.17).).  The  reference  spectrum  solution 
is  also  obtained  by  using  the  boundary  conditions  (2.18) 
to  (2.2l).  For  r  <  £  *  the  boundary  conditions  on 
the  solution  of  equation  (4.13)  are 

u g  (c)  =  0  (4.14) 


and 


dr  (u7r)  •  rJi(kr)  ) 

d  ,  hb>f1  .  4 

dr  (i7r)  ) 

Ug(r)  -  rj^(kr) 

r=£ 

rh^1)  (i7r) 

(4.15) 
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The  separation  distance  is  defined  by  equation  (4.9). 

This  condition  may  be  expressed  in  partial  waves  by  an 
equation  similar  to  (2.23)  except  that  the  statistical 
factors  c ^  must  be  used.  Thus 

£(k) 

4tt  2  (2 i  +  l)  c %  J  r  j^(kr)  v(r)  u^(r)  dr  =  0  (4.l6) 


We  can  substitute  in  the  integrand  for  v(r)u^(r) 
from  equation  (4.13),  then  use  integration  by  parts  and 
substitute  boundary  condition  (4.15).  The  result  of 
this  manipulation  is 


0 


47 T 


* 


m 


2  (2 £  +  1)  c  n 
l  1 


r 

< 


i 

(y2  +  k2)  /  rj^(kr) 

0 


[rj^ (kr)-u^ (r) ]  dr 


+ 


u/g(^)-r’J^(kr) 

r  h^1 ^  (iyr) 


(  r j ^ ( kr )  rh^1 ^ (iyr  )- 


rh^1  ^  (iyr)  r^(kr)} 


(4.17) 


The  separation  distance  |R  for  the  reference  spectrum 
approximation  is  calculated  by  using  the  first  order 
approximation  £>q  to  obtain  the  first  approximation 
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to  the  energy  denominator  e^k,  k*)«  Then  eR  (k,k») 
is  obtained  by  fitting  e^  with  the  quadratic  form 
(4.10)  (the  best  convergence  is  obtained  by  choosing  the 

*  2  /  x 

parameters  m  and  y  (k)  to  approximate  as  closely 
as  possible  the  actual  energy  spectrum  for  k*  between 
2  kp  and  4  k^  according  to  Bethe  et  al^11^). 

Equation  (4.1l)  is  then  solved  for  ip  subject  to  the 
given  boundary  conditions  and  a  second  order  approximation 

ID 

to  £  obtained  from  equation  (4.17).  The  process  is 
repeated  until  self-consistency  is  achieved. 

It  is  now  necessary  to  discuss  corrections 
to  the  reference  spectrum  result.  We  put 

|  =  £R  +  A  Ix  +  &§p  +  A?s  (4.18) 

where  the  correction  terms  are  the  Pauli  correction, 
center  of  mass  correction  and  spectral  correction, 
respectively.  These  terms  correct  the  inaccuracies 
due  to  the  three  reference  spectrum  approximations 
listed  previously.  The  reference  spectrum  approximation 
to  the  GS  matrix  is  easily  obtained  from  the  solutions 
£R(k)  and  ipR{tf,r) .  The  result  is  given  by 
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(4.19) 


and  the  equation  satisfied  by  the  GRs  matrix  is 


GRs  =  vRs  +  vRs  1  gRs 

n 


(4.20) 


e 


The  GS  matrix  itself  satisfies 


(4.21) 


where  Q  and  e  depend  on  the  parameters  and  P. 

r 


If  only  the  Pauli  correction  is  to  be  made 


the  energy  denominators  will  be  the  same  in  equations 
(4.20)  and  (4.2l).  The  correction  for  a  non-zero 
center  of  mass  momentum  will  involve  a  change  in  the 
energy  denominator  and  in  the  exclusion  operator  Q. 
The  effect  is  most  important  in  the  exclusion  operator 


Q(it?)  =  e(q2  -  kp2)  e(kj2  -  kp2) 


where  a  =  cos  (Ic, "?).  It  is  convenient  to  average  over 
the  angle  by  putting 
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Q(k,P)  =  Q(k,P) 


\  J  da 

-1 

o 

2  P  2 

o  if  c  kpd 


k2  + 


2  2  P  2 

if  kp  <kN-^-  < kP+k^ 


=  1 


if  kP  +  k. 


2  .  .  2 
i  ^  k  + 


P2 

w 


(4.23) 


The  change  in  the  energy  denominator  is  small  and  is 
neglected.  Now  on  comparing  the  two  equations  (4.20) 
and  (4.21)  with  e  =  e  we  have 


G  =  G 


Rs  -  vRs(l  +  i  GRs)  +  Vs  (1  +  S  Gs) 


,Rs 


=  G  -  v 


Rs  (1  +  1  qRs )  (1  +  a  Qs)  +  GRs  S  GS 
'  e  e  e 


Q  ns  ( i  _i_  iL  nRs  _  os  2.  qRs 


+  vb  (l  +  —  Gb )  (1  +  ±  Gnb)  -  Gk 


nRs  s  Rs  ,  /  s  Rs \  / 1  nRs  Q  rs \  ,  nRs 

G  +  v  -  v  +  (v  -  v  )  {—  G  +  —  G  J  +  G 

e  “ 


-  GS  i  GRs  +  vs  S  GS  1  GRs  -  vRs  i  GRs  $  Gs 

e  e  e  e  e 


A  first  order  solution  is  obtained  by  neglecting  the 
terms  of  third  order  in  v  and  G  and  by  approximating 


G?|(D 
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s  s 

v  and  G  by  the  reference  spectrum  results  in  the 
second  order  terms.  We  are  then  left  with 


GS  2  G 


Rs 


+  vk 


vEs+  qRs  Q-i  gRs 


(4.24) 


c*  p  o 

The  diagonal  elements  of  both  G  and  G  vanish 

p 

according  to  the  definitions  of  £  and  |  so  that 
for  k  < 


<i?  I  v 


s  Rs 

-  v 


i  k  >  = 


(277-) 


<C  I  GRs  |  £'>[l-Q(k',P,k  )] 
1  /  dl  -  F 


3 


eF  (k,k') 


<£*  |  GRs  |  ic  > 


(4.25) 


Rs 

The  matrix  elements  of  G  are  given  by  equation  (4.19) 
so  the  right  hand  side  of  equation  (4.25)  is  a  known 
function  of  k,  P  and  which  shall  be  denoted  by 

P  (k,  F,  k  )  for  convenience.  Then 


F(k,P,kp)  =  f  e  lk#rv(r) [@(^-r)-0(^R-r)]  elk  r  dr 

£  2 
=  4t t  J  v(r)  r  dr 

<=7  4tt  v(£R(k)  )  UR(k)  ]  [5(k^P,kF)-?R(k)  ] 


(4.26) 
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The  solution  [£(k, P,kp)  -  |R(k)]  is  the  sum  of 
the  correction  terms  A£  and  A^p  of  equation  (4.18). 

The  correction  terms  could  be  handled  separately,  but 
the  correction  due  to  the  exclusion  principle  alone, 

A£  .>  is  just  the  solution  of  equation  (4.26)  for  the 

u\. 

special  case  P  =  0. 

The  spectral  correction  is  known  to  be  small 
A  first  order  correction  would  be  calculated  by  comparing 
equation  (4.20)  and  (4.2l)  with  Q  set  equal  to  unity. 
The  equation 


Rs  .  ~Rs 
v  +  G 


,1  lx  nRs 

'  *~R '  G 
e 


(4.27) 


is  obtained  by  a  derivation  similar  to  that  used  for 
equation  (4.24).  Then  on  taking  a  diagonal  element  we 
have 

<k  |  vS-vRs  J  dkT  <k  |GRs  |  kT>(  R- - - —  - 

(2t r)3  e  (k,kT) 


'  e(k <k’  I  °RS  I  k  >  (4‘28) 

For  e (k, k1 )  we  use  the  expression  obtained  by  putting  £  (k) 
in  equations  (3-40) ,  ( 3 • 51 )  and  (3.4l).  The  spectral 
correction  is  then  calculated  from  equation  (4.28)  by  the 
same  procedure  as  was  used  in  obtaining  corrections  for  the 
exclusion  principle  and  center  of  mass  motion  from  equation  (4.2b). 
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Chapter  5.  Divergence  of  the  Gs  -  matrix  Expansion 


It  is  well  known  that  the  Brueckner  reaction 
matrix  expansion  is  divergent  for  matrix  elements  taken 
sufficiently  close  to  the  fermi  surface.  In  this  chapter  it 
is  shown  that  these  difficulties  are  also  present  in  the 
expansion  of  the  G  -  matrix  defined  earlier.  There  will 
follow  a  discussion  of  efforts  to  modify  the  theory  so 
as  to  remove  this  difficulty. 


The  integral  equation  for  G  is 


ns  s  .  s  Q, 

G  =  v  +  v  —  G 


(5.1) 


Restricting  the  discussion  to  the  case  where  the 
center  of  mass  momentum  is  zero  we  have 


|  G  |  Icq> 


=  <Tc  |  vs  | 


_y  1  , 

k  >  +  - 5-  /  dk  — 

°  l»V>K„ 


Vs  I  I?  ><£■'  |GS  |  to  > 


e (k  , k' ) 


o 


(5.2) 


where  It  it  since  Gs  is  only  defined  by  its  action 

on  the  ground  state  wave  vector.  It  is  now  assumed 

that  the  potential  can  be  expanded  in  spherical  harmonics 

(20) 


in  the  form 
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im(?)  W?,)  (5-3) 

1=0  m =-£ 

The  local  potential  used  previously  can  be  obtained  by 

O 

using  v^(r,rT)  =  v(r)  8(r-r,)/r  .  At  this  point 

we  want  to  illustrate  the  divergence  of  the  expansion 

g 

of  G  and  to  this  end  we  assume  a  less  realistic 

( 21 ) 

separable  potential v  ’  by  putting 

vi  ( r  ■>  r  1  )  =  (r)  f^(rT)  (5.4) 


<r  |v  fr  >  = 


vhr'r 


)Y 


The  short  range  potential  and  separation  distance  are 
now  introduced  by  putting 


v 


S'  -  '  * 


i 


(r,r’  )  =  «f^  (r)  f^r1)  Q(i-r') 


l 


* 


=  otf^  (r)  f^(r’)  . 


(5-5) 


Then 


00  i 


tvs  I  ?  >  =  2  2  af*(r  )f*(r  )  Y^(?)  *Jm(*'),(5.6) 

£=0  m=-£ 

and  the  matrix  elements  <k  |  v  |  kl>  are  given  by 

<k  |  vS  |  ic*>  =  / /dr  dr1  (ic,  r)  <r  |  vS  |  r  1  >  <!>  (k1 ,  r 1 ) 

=  2  2a/  dr  4>*(ic,r)  ffr)Yjm(r)  /dr4>  (T?> ,  r  )f|  (r)  Y^m(r) 

l  m 


(5.7) 
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It  is  now  assumed  that  <f>  (5T,  r)  can  be  approximated  by 
a  plane  wave  and  a  partial  wave  decomposition  is  used 
to  write 


*(k,  ?)  =  e1  k'r  =  4tt  2  Ei1  j^(kr)  Y^m(k)  Y^m(?) 

1  m 


(5.8) 


Then  with  the  orthonormality  property  of  the  spherical 
harmonics  it  follows  that 


00 


/dr  0(k,r)  f^(r)  Y^(r)  = 


4tt  i^  /  j^(kr)  f^(r)  r2dr  Y^(k) 


0 


im 


and  equation  (7)  becomes 


* 


<k|vs|k’>=  2  ahjg(k)  h^  (k T )  Y^m(k)  Y^Jk') 


i,  m 


(5.9) 


where 


and 


00 


h^(k) 


=  4t t  f  j^(kr)  fg(r)  r  dr 


0 


00 


h»  (k)  =  4t tJ  j„  (kr)  f®(r)  r2dr 


7 


0 


1,  v  '  \e 


(5.10) 


=  4tt  /  j^(kr)  f^(r)  r  dr 


0 


(5.11) 


g 

For  the  G  matrix  we  also  assume  an 
expansion  in  spherical  harmonics  of  the  form 


<lC  I  GS  |  k_> 


=  2 
&  9  m 


(k,k  J  Y*q  (k) 
'o'  Jom  ' 


Y 


im 


(k  ) 


o 


(5.12) 


Equations  (5*9)  and  (5.12)  may  now  be  substituted  in 
equation  (5.2).  The  following  integral  equation  for 

g 

G^(k,k  )  is  obtained 


G?k>ko> 


*,  .  S/  .  “  Mk)  “  h|  k>)Gf  k',k 

a  h«k)  h®  k  +  - /  — - * - °  k,2dk. 

1  °  (2^)3  k  e(k  ,k>) 

o 


(5.13) 


This  can  be  solved 


exactly 


(22) 


The  result  is 


G^(k,ko) 

where 


a  h](k)  (kQ) 
1  -  “  Fi  ^ 


I 

kF 


h](k)  h® (k) 
e (kQ,k) 


k^dk. 


(5.14) 


(5.15) 


Since  the  separation  distance  does  not  exist  for 
k  >  kp  ,  h^(k)  ls  eGual  to  h^(k)  in  the  reSlon  of 
integration  in  equation  (5.15).  Therefore 


= 
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(kj  =  F„(kJ  = 


o 


i  v  o 


(27 r) 


3 


00 

/ 

k 


h,(k)  |  2 

^  k  dk 


F 


e(kQ,k) 


(5.16) 


Bell  and  Squires  have  shown  the  divergence 
of  the  reaction  matrix  for  a  separable  potential  by 
a  method  which  is  adopted  here  for  the  Gs  -  matrix. 

It  is  assumed  that 


(1)  h^(k)  does  not  vanish  in  the  region  k^<k<oo  and 

(2)  e(k  ,k)  is  negative  for  all  k^k^  k>kp. 


Then  the  integrand  in  equation  (5.1 6)  is  always  negative 
and  for  any  positive  A 


<  ,dr  / 


kF+  X 


h£(k) 


(2t r) 


k 


F 


e(kQ,k) 


k^dk 


(5.17) 


For  any  reasonable  form  of  e(kQ^k)  there  exists 
some  negative  number  A  such  that 


1  hbk)  1  <  a 

e ( kQ ,  k )  k  -  kQ 


(5.18) 


for  k  <k„,  k„<k<k„+A.  Then  it  follows  that 

O  r  r  r 


A  dk 


=  A  log 


k^  +  A  -  k 
F  o 


Vko}  < 


kp+A 

/ 

kF 


k-k 

o 
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(5.19) 


Thus  it  is  proved  that  F  (k  )  ->  -  »  as  k  k„. 

Jo  o  O  -t1 

The  series  expansion  of  Gs  may  be  written  in  terms  of 

F„(k  )  as  follows 
i'  o 


<£  I  g! 


I  Ic  > 
1  o 


=  2 

&  9  m 


* 


a  h^(k) 


Vko} 


*  . 
Wk> 


Y.  (k  ) 
&  m  o 


[l  +  a  F 


i 


(k  ) 
v  o  ’ 


+ 


a2F 


i 


(k  ) 
v  o ' 


+ 


(5.20) 

.  ] 


Therefore  no  matter  what  the  sign  and  magnitude  of  a, 

g 

the  series  expansion  of  G  diverges  for  a  range  of 

values  of  k  near  the  fermi  surface. 

o 


In  addition  it  should  be  noted  that  for 

negative  a,,  the  denominator  of  equation  (5.14)  can 

vanish  unless  the  minimum  values  of  a  F„(k  )  exceeds 

unity  (if  this  is  true,  expansion  (5.20)  diverges  for 

all  k  ).  In  general  a  Ffl(k  )  will  vanish  for  some 
o'  0  i v  o ' 

k  and  at  this  value,  the  reaction  matrix  will  have 
o 

a  singularity.  For  any  given  i,  the  singularity 
occurs  at  the  lower  bound  of  the  region  of  divergence 
of  the  contribution  to  equation  (5.20)  of  the  i ^  partial 
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wave.  .For  a  finite  system  the  singularity  may  not  occur 
for  the  allowed  discrete  values  of  k  ,  depending  on 
the  volume  of  the  system  and  the  strength  a  of  the 

g 

interaction.  As  these  are  varied,  G  oscillates 
infinitely. 


The  existence  of  singularities  in  the  reaction 

matrix  and  the  divergence  of  the  reaction  matrix 

expansion  are  not  phenomena  peculiar  to  separable  potentials. 
(23 ) 

Emery x  '  obtained  similar  results  for  a  local  potential. 

His  method  is  adopted  here  for  the  investigation  of 
the  divergence  of  the  Gs  -  matrix.  This  time  it  is 
more  convenient  to  work  with  the  integral  equation  (3.2 8) 
for  the  wave  function.  With  the  assumption  that  the 
center  of  mass  of  the  two  interacting  nucleons  is  at 
rest,  the  two  particle  wave  function  in  the  co-ordinate 
representation  satisfies  the  equation 


^(k,r)  =  $(Ie,r)  +  /  K  (r,  r 1  )  vs  (r*  )  ^(It,r’T)  dr”1  (5*21) 


where 


_ 1_ 

(  2tt) 


3 


$  (ic» ,  r)  <t>  (it1 ,  r 1 ) 

I  • 

kT>kp 


die 1 


K(r’,r’1 )  = 


e (k, kr ) 


(5.22) 
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To  simplify  the  equations  we  assume  that  v(r)  is  a 
central  potential  and  we  use  the  plane  wave  approximation 
for  the  unperturbed  wave  function  4>  and  then  use 
the  S  -  state  approximation ,  i.e., 

.  ->n  sin  kr 

0  (k’r)  =  kr 

Furthermore  only  kinetic  energy  terms  are  included 
in  the  energy  denominator,  i.e.. 


e  (k,  k T )  =  k^  -  k  *  ^ 


With  these  crude  assumptions  equations  (5.21)  and 
(5.22)  become 


00 


r  f  (k,  r ) 
where 


sin  kr 
k 


+  /  K(r,rT)  vs(r!)  ^(k,r»)  rTdr*  (5.23) 

0 


00 


K  (r, r ’  )  = 


2  r  sin  kTr  sin  k’r1 
k^ 


^  / 
IT  J 


2  2 

k  -  kl 


dk* 


The  series  solution  for  r f  (k,r)  is  then 


(5.24) 


00 


r  f  (k,  r )  = 


sin  kr  .  r  \  _S/_.\  sin  krT 


k 


+  f  drT  K(r,r!)  v  (rT) 


0 


k 


+ 


(5.25) 


This  is  a  Neumann  series  which  gives  the  solution 
provided  that  the  series  converges  uniformly.  A  theorem 
of  Fredholm  (see,  for  example,  reference  (24)  )  states 
that  if  =  1  is  an  eigenvalue  of  the  corresponding 
homogeneous  equation 


r  ^(k,r)  =  Ai(k,kF)  /  K(r,rr)  vs(rT)  rl  ip±(\zJrt)  dr1  (5.26) 

0 

then  the  nonhomogeneous  integral  equation,  (5*23)* 
has  a  solution  r  ^(k,r)  if  and  only  if 


00  •  i 

r  .  /,  s  sin  kr  ,  „ 

J  r  ^(k,r)  — ^ -  dr  =  0 

0 

But  it  was  earlier  specified  that  the  normalized 
many-body  wave  function  |  >  and  the  unperturbed  wave 

function  I  0  >  satisfy  the  relation  <  0  |  I  >  =  1 . 
Therefore  the  corresponding  two-particle  wave  functions 
can  not  be  orthogonal.  If  A^(k,kp)  =1  is  an  eigenvalue 
of  equation  (5.26),  then  equation  (5.23)  has  no  solution 
and  the  Neumann  series  (5.25)  is  divergent.  Furthermore 
if  A  (k,kp)  is  any  eigenvalue  of  equation  (5.26), 
the  equation 


sin  kr 


+  ~K 


f  K(r,rT)  vs(rT)  r*^(k,rT)drT 
0 


r  if/  (k,  r )  = 


k 
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has  no  solution  and  the  associated  expansion, 

00 

r  ^(k,r)  =  Sl^  kr  +  A  /  K(r,r«)  vs(r»)  S1^.-kr.t..  dr* 

0 

0°  00  *t 

+  A  /  dr1  /  dr**  K(r,r*)  vs(r*)  K(  r * ,  r"  -  + 

0  0 

is  divergent.  Then  if  |  A  |  <  1  ,  the  expansion  (5.25) 
must  also  be  divergent.  Thus  if  there  are  eigenvalues 
of  equation  (5.26)  satisfying  |  Ai(k,kp)  1, 
the  series  of  equation  (5.25)  will  diverge  and  if 
A^(k,kp)  =  1  is  an  eigenvalue,  the  integral  equation 
is  singular. 

g 

The  series  expansion  of  G  is  obtained  by 
substituting  the  expansion  of  if/(\z,r)  into 

<Ic  |  Gs  |  l?o>  =  /  dr  <t>*(£,r)  vs  ( r )  r)  , 

g 

so  that  the  expansion  of  G  has  the  same  properties  as 
the  expansion  of  the  wave  function  f  (Ic,?).  Thus  the 

g 

problem  of  the  existence  of  singularities  in  G  and 

the  divergence  of  the  associated  series  expansion  is 

reduced  to  a  discussion  of  the  eigenvalues  of  equation  (5.26). 


The  following  properties  of  the  eigenvalues 
V  (k,kp)  were  derived  by  Emery '  D J  on  the  assumption 
that  the  matrix  elements  of  the  potential  <k  |  vs  |  kT> 
are  finite  (an  infinite  hard  core  could  then  be 
considered  as  the  limiting  case  of  a  finite  core).  The 
derivation  will  not  be  repeated  here. 


I. 


II. 


III. 


IV. 


V. 


The  eigenvalues  are  real  and  V (k,kp)  i=  0  unless 

k  =  k^. 

d  I  K  (k,k  )  | 

-  >  0  for  0  <  k  ^  k„  and  all  i. 

Up  F 

d  I  V  (k,k  )  | 

-  ^  0  for  0  ^  k  <  k  and  all  i. 

d  k  ^ 

a)  If  <kp  |  vs  |  kp>  <  0,  An  -o‘  as  k-*kp 
for  one  and  only  one  value  n  of  i; 

b)  If  <kp  |  vs  |  kp>  >  0,  An  -0"  as  k-kp 
for  one  and  only  one  value  n  of  i; 

c)  If  <k„  |  vs  |  k  >  =  0,  0  as  k-4c  for  at 

least  2  values  of  i.  In  general  there  are  two 
such  values,  one  tending  to  0+,  the  other  to  0  . 

a)  If  vs(r)  has  no  bound  states,  there  is  no 
eigenvalue  at  k=kp=0  which  satisfies  0<V  (0,0)<1; 

b)  If  vs(r)  has  N  bound  states  there  are  N  values  of 
i  for  which  0<V  (0,0)^1  and  equality  holds  only  if 


there  is  a  bound  state  of  zero  energy. 
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That  these  properties  of  the  A  imply 

i 

divergence  of  the  series  expansion  of  G  may  be  seen  as 

follows.  For  any  e  >  0,  it  follows  from  III  and  IV 

that  there  exists  a  value  k  <  k^  such  that 

o  F 

I  \(k,kp)  |  <  € 

for  at  least  one  value  of  i  and  for  all  k  in  the 
interval  k  ^  k  .  Thus  there  is  always  a  finite 

O  r 

range  of  momenta  (with  kp  as  the  upper  bound)  for  which 
at  least  one  eigenvalue  of  Eq.  (2 6)  has  absolute  value 
less  than  one  and  for  which  the  series  expansion  of 

g 

G  therefore  diverges.  This  is  identical  to  the  conclusion 
reached  for  the  case  of  a  separable  potential. 

We  also  want  to  examine  the  conditions  for  the 

g 

existence  of  a  solution  to  the  integral  equation  for  G  . 

We  know  that  the  equation  is  singular  if  some 
A. (k,k„)  =  1.  From  property  III  of  the  eigenvalues  we 

1  r 

seek  a  A^  such  that  A^(0,kp)  >  1  while  A^(kp,kp)  <  1. 

Q1 

Now  since  v  (r)  has  no  bound  states.,  it  follows  from 
V  and  II  that  all  positive  eigenvalues  satisfy 
A^(0,kp)  >  1.  Then  assuming  the  existence  of  positive 
eigenvalues,  a  singularity  will  exist  if  A1(kpkp)  <  1 
for  one  of  these  positive  eigenvalues.  Now  if  the 
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potential  vS(r)  is  replaced  by  \i  v3 (v) ,  the  matrix 

g 

G  will  be  singular  if  one  of  eigenvalues  A.(k,k  ) 

1  r 

is  equal  to  \±.  Then  the  necessary  condition  for  a 
singularity  arising  from  A‘.  is  p  >  A.  (k^k^)  where 

1  1  r  r 

is  a  positive  eigenvalue.,  i.e.  the  potential  must 
be  sufficiently  strong. 


There  is  one  eigenvalue  A  for  which 

n 

An(kp;,kF)  =  0.  Assuming  <kp  |  vs  |  kp>  <:  0  (Property  IV) 

this  eigenvalue  can  lead  to  a  singularity  with  a  potential 

[i  vs(r)  only  if  [i  <  A  (0,k  );  i.e.  if  the  potential  is 

n  F 

sufficiently  weak.  This  type  of  singularity  is  called 
a  Cooper  singularity.  There  is  no  Cooper  singularity  if 
<kF  |  vs  |  kF>  >  0  since  A^(k,kF)  of  property  IV  is 
then  negative  everywhere.  This  is  the  case  in  nuclear 
matter  where  an  arbitrarily  strong  repulsive  core 
exists . 


( 2A ) 

Emery v  J  solved  the  eigenvalue  equation  (5.26) 

( 2R ) 

using  the  Gammel-Thaler  potential v  ^ '  for  singlet  and 
triplet  states.  The  lowest  eigenvalues  A^k^kp)  were 
found  to  be  0.56  (singlet)  and  0.75  (triplet).,  both 
eigenvalues  leading  to  singularities.  The  result  is 
not  very  meaningful  in  view  of  the  crude  approximation 
of  neglecting  the  potential  energy  terms  in  the  energy 
denominator  of  equation  (22).  A  more  common  approach  is 
to  use  the  effective  mass  approximation  for  e(kjk’ ), 
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i.e.  to  put 

e(k,k*) 


* 

m 


With  this  approximation  the  Gs  matrix  is  singular  if 

some  A. (k^k^j  =  m  .  Hence  a  value  of  m*  less  than 
1  F 

0.56  or  0.75  would  be  necessary  to  remove  the  singlet 

or  triplet  singularities  respectively.  In  the 

calculation  of  the  reaction  matrix  by  Brueckner  and 

Gammel^  ,  which  used  the  Gammel  -  Thaler  potential, 

■x- 

a  self-consistent  potential  was  obtained  with  rn  =  0.73 
at  the  Fermi  surface  and  increasing  with  k  so  that  the 
singlet  singularity  is  not  removed  and  the  triplet  singu¬ 
larity  would  most  probably  not  be  removed.  The  singularities 
were  not  noticed  in  the  calculations  of  Brueckner  and 
Gammel  because  of  inaccuracies  in  their  numerical 
approximations.  In  general  a  numerical  solution  of  the 
integral  equation  is  unlikely  to  reveal  this  effect 
because  it  is  important  only  for  k  very  close  to  the 
Fermi  surface. 

g 

The  divergence  of  the  G  matrix  for  values  of 
k  near  k^  is  due  to  the  fact  that  for  states  with 
momenta  sufficiently  near  the  Fermi  momentum  the  energy 
denominator  becomes  arbitrarily  small  (in  the  effective 
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mass  approximation)  compared  with  the  strength  of  the 
interaction.  A  different  form  for  the  single  particle 
potential  U(k^)  would  not  be  expected  to  eliminate  the 

g 

divergence  of  G  for  large  k  in  the  range  o  «  k  «  k 
as  long  as  U(k^)  is  continuous  at  the  Fermi  surface 
because  then  the  energy  denominator  could  still  be 
arbitrarily  small.  One  can  however  remedy  the 
situation  by  assuming  that  U(k^)  is  a  monotonic 
increasing  function  of  with  a  discontinuity  at  the 
Fermi  surface.  The  possibility  of  an  energy  gap  suggests 
that  a  "superfluid"  state  exists  for  nuclear  matter. 

One  could  then  hope  to  apply  the  Bardeen  -  Cooper  - 
Schrieffer  (BCS)  theory^  '  to  the  problem  and  obtain 
an  energy  gap  in  a  logical  fashion.  However  it  is  not 
by  any  means  clear  that  the  use  of  the  BCS  model  is 
justified  for  nuclear  matter.  The  hard  core  of  the 
nucleon  -  nucleon  interaction  would  demand  modification 
of  the  BCS  trial  wave  function  and  the  BCS  method  does  not 
contain  any  provision  for  dealing  with  the  self-consistency 
requirements  encountered  in  the  study  of  nuclear 
matter.  Attempts  to  demonstrate  the  superfluidity  of 
nuclear  matter  have  so  far  been  inconclusive  (see  review 
in  Bell  and  Squires'').  Those  who  have  succeeded  in 
predicting  an  energy  gap  on  the  basis  of  the  BCS  theory 3 
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have  obtained  gaps  of  about  0.1  Mev  or  less  which 

would  not  be  large  enough  to  eliminate  the  singularity 

( 2Q ) 

and  the  region  of  divergence.  Nuttall  et  al  have 

shown  that  with  the  Gammel  -  Thaler  potential.,  the 
singularity  in  the  Brueckner  reaction  matrix  disappeared 
when  the  gap  was  about  0.4  Mev  or  more. 

The  inclusion  of  additional  terms  in  the 

perturbation  series  might  be  expected  to  lead  to  a 

non-singular  equation  for  the  reaction  matrix.  Recall 

that  the  Hamiltonian  defined  earlier  included  only 

particle-particle  interactions.  The  series  expansion 

thus  contains  only  so-called  ladder  diagrams  and  is 

referred  to  as  the  "ladder  approximation".  A  more 

exact  theory  would  also  include  the  effects  of  hole-hole 

interactions  and  "pair  production"  interactions.  However 

it  was  shown  by  Mehta ^0)  that  with  a  separable  potential 

the  reaction  matrix  defined  by  taking  hole-hole  interactions 

as  well  as  particle-particle  interactions  has  the  same 

singularity  as  the  Brueckner  reaction  matrix.  Bell  and 

( 7) 

Squires  v  /  included  the  pair-production  interactions  as 
well  and  came  to  the  same  conclusion.  This  was  also  based 
on  the  separable  potential  model.  Thus  it  would  seem  that 
the  ladder  approximation  is  not  the  source  of  the 
singularity . 
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Chapter  6.  Conclusions 

The  preceding  chapters  have  outlined  a 
method  of  evaluating  the  Brueckner  reaction  matrix. 

This  method  is  an  improvement  over  a  direct  solution 
of  the  Brueckner  equations  because 

a)  the  short  range  two-nucleon  interaction,  which  is 
not  well  understood,  is  treated  as  a  perturbation 

and  is  not  required  for  a  first  order  calculation,  and 

b)  the  self-consistency  requirements  on  the  final 
solution  are  satisfied  in  the  calculation  of  the 
separation  distance  where  the  mathematics  is  simplified 
by  the  use  of  the  reference  spectrum  approximation. 

This  approach  to  the  problem  is  quite  similar 

( 9) 

to  the  separation  method  of  Moszkowski  and  Scott v  ' . 
However,  as  mentioned  earlier,  the  Moszkowski  -  Scott 
separation  is  so  defined  that  the  short  range  interaction 
gives  zero  phase  shift  for  free  particle  scattering  and 
there  will  be  a  small  first  order  contribution  to  the 
energy  from  the  short  range  interaction.  A  more  important 
criticism  of  the  Moszkowski  -  Scott  method  is  that 
their  separation  distance  depends  on  the  angular 
momentum  H  of  the  interacting  pair  as  well  as  the 
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relative  momentum  k.  Thus  in  states  such  as  and 

3 

P-^  where  the  potential  is  on  the  average  repulsive  the 

( ) 

separation  distance  does  not  exist v  .  With  our  definition 
the  existence  of  the  separation  distance  is  guaranteed 
for  relative  momenta  less  than  or  equal  to  the  fermi 
momentum  if  the  ground  state  of  the  system  is  bound. 

In  addition  the  numerical  results  of 
Moszkowski  and  Scott  show  that  the  direct  contribution 
to  the  reaction  matrix  from  the  short  range  interaction 
is  positive  and  cancels  part  of  the  long  range  contri¬ 
bution.  Hence  with  our  formulation  a  larger  separation 
distance  is  expected.  The  long  range  potential  is 
thus  weaker  and  we  would  expect  a  better  convergence 
of  the  expansion  of  G  and  a  more  accurate  first  order 
result . 

Bethe  et  al^^  have  discussed  a  separation 
within  the  framework  of  the  reference  spectrum  method. 

They  separate  the  interaction  at  the  point  where  the 
reference  spectrum  wave  function  joins  smoothly  with 
the  unperturbed  wave  function.  This  definition 
leads  to  a  very  large  direct  contribution  to  the  energy 
from  the  short  range  interaction.  In  general  it  can 
be  said  that  no  matter  how  the  separation  distance  is 
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defined,  the  plane  wave  approximation  leads  to  an 
expansion  given  by  equation  (3.34).  It  is  reasonable 

Cj 

to  suppose  that  the  off-diagonal  elements  of  G  are  on 
the  whole  smallest  when  the  diagonal  elements  vanish 
and  thus  the  best  convergence  should  be  obtained  with 
our  definition  of  the  separation  distance. 

Another  advantage  of  this  formulation  of  the 

problem  is  the  possibility  of  using  the  Hartree  -  Fock 

technique  to  evaluate  the  eigenfunctions  of  the  effective 

Hamiltonian.  This  makes  the  method  useful  even  if  the 

tensor  force  is  strong.  The  Hartree  -  Fock  method  is 

(l4 ) 

quite  complicated  and  details  are  given  elsewhere v  . 

The  calculations  are  much  simpler  if  one  uses  the  plane- 
wave  approximation. 


The  singular  nature  of  the  integral  equations 

of  the  Brueckner  theory  has  not  been  mentioned  in  chapters 

3  and  4.  This  point  is  discussed  in  chapter  5  and  the 

possibility  of  a  gap  in  the  single  particle  energy 

spectrum  is  introduced.  This  would  require  modification 

of  the  energy  denominator  by  the  introduction  of  the 

"off  the  energy  shell"  definition  of  the  single  particle 

(8) 


potential  energy  for  excited  states 

.  (12) 
shown '  ' 


It  can  easily  be 
that  an  energy  gap  is  readily  incorporated  in 


the  reference  spectrum  method. 
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No  numerical  calculations  have  been  included 

to  test  the  method.  It  would  be  of  interest  to  compare 

numerical  results  with  those  of  Brueckner  and  Gammel^^ 

(g) 

and  Moszkowski  and  Scott'  ' .  However  it  is  of  little 
consequence  whether  or  not  we  get  results  in  closer 
agreement  with  experiment.  Comparison  with  experiment 
shows  whether  the  assumptions  made  about  nuclear 
interactions  are  acceptable  or  not.  As  Bethe^1^  puts 
itj  "The  method  must  rest  entirely  on  mathematical 
proof  and  internal  consistency. " 
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